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1. Introduction 

Supersymmetric field theories play a central role in modern theories of particle physics. 
They offer a possible solution to the gauge hierarchy problem and are often more tractable 
analytically than their non-supersymmetric counterparts while still exhibiting features like 
confinement and chiral symmetry breaking [||. Super Yang- Mills theories are especially 
interesting because of their possible connection to string and M-theory ||. One example 
of this is the conjectured equivalence between (p + l)-dimensional super Yang- Mills theory 
and supergravity containing Dp-br&ne sources 

One attractive scenario for embedding the Standard Model in such a theory is to 
imagine that supersymmetry breaks spontaneously as some gauge coupling becomes large at 
low energy. Unfortunately various non-renormalization theorems forbid such a spontaneous 
breaking at any finite order of perturbation theory so one must turn to non-perturbative 
mechanisms to drive such symmetry breaking. 

Of course the lattice furnishes perhaps the only generally applicable way to study non- 
perturbative dynamics in field theory. However, the difficulties of discretizing supersym- 
metric theories are well known. Generic naive discretizations of continuum supersymmetric 
theories do not preserve supersymmetry. Typically quantum corrections then generate a 
large number of relevant supersymmetry violating interactions whose couplings must be 
tuned to zero as the lattice spacing is reduced. This is both unnatural and in many cases 
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(especially for models with extended supersymmetry) prohibitively difficult. Various at- 
tempts have been made over the last twenty five years to overcome these problems see § 
and H and the recent reviews jj], ||, ||]. 

Quite recently a series of new approaches have been developed which share the common 
feature of preserving a sub-algebra of the full supersymmetry algebra exactly at finite lattice 
spacing 1 In the approach pioneered by Kaplan and collaborators [pd| the lattice theory is 
derived by applying a carefully chosen orbifold condition to a supersymmetric matrix model 
(see also []l2|, [L4|, |l5| ) . A second approach proceeds from a reformulation of the theory 
using ideas drawn from topological field theory [16]. . It is this approach that is the focus 
of this paper. The hope is that the residual exact supersymmetry will protect the theory 
from dangerous radiative corrections and obviate the need for fine tuning [11, [l8|, 17 1. 



This technique was initially used for theories without gauge symmetry [T^, 2C, ^] cor- 
responding to supersymmetric quantum mechanics, the 2D complex Wess-Zumino model 
and supersymmetric sigma models. An implementation for gauge theories was initially 



given by Sugino [22|. Important progress was made when Kawamoto and collaborators 
|2^ | pointed out the connection between these topological formulations and Kahler-Dirac 
fermions. This has led to additional lattice formulations which emphasize the geometrical 
nature of the underlying theory [[h]] and [l], |24j] (Hamiltonian formulations of lattice super- 



symmetric theories employing Kahler-Dirac fermions were first proposed in [25| and later 
in H) 

The key starting point in this approach is to construct a new rotation group from a 
combination of the original rotation group and part of the R-symmetry associated with the 
extended SUSY. The supersymmetric field theory is then reformulated in terms of fields 
which transform as integer spin representations of this new rotation group 



process is given the name twisting and in flat space one can think of it as merely an exotic 
change of variables in the theory. In this process a scalar anticommuting field is always 
produced associated with a nilpotent supercharge Q. Furthermore, as argued in |^8L [l]] the 
twisted superalgebra implies that the action rewritten in terms of these twisted fields is 
generically Q-exact. In this case it is straightforward to construct a lattice action which is 
Q-invariant provided only that we preserve the nilpotency of Q under discretization. 

In (l[ we gave a concrete application of these ideas by constructing a gauge invariant 
lattice regularization of M = 2 super Yang-Mills theory in two dimensions which was invari- 
ant under the twisted scalar supersymmetry and an additional SO(l, 1) scaling symmetry. 
In this paper we present results from the first numerical simulations of this model. 

In this initial work we have examined a variety of Ward identities following from the 
twisted supersymmetry. In the case of U (2) we find generally good agreement with theory 
for sufficiently large coupling j3. Specifically we do not need to fine tune any additional 
couplings to see a restoration of the twisted supersymmetry for small lattice spacing. In 
the case of SU{2) the agreement is even better - we see no statistically significant breaking 



x Very recently a lattice construction of Af = 2 super Yang-Mills in D = 2 has been proposed which 
claims to preserve all the supercharges [[To| 

2 It appears that these two approaches may be intimately connected - private communication Mithat 
Unsal 



- 2 - 



of the Ward identities even for small (3. Our results for the string tension indicate that the 
lattice theory possesses a single confining phase. In addition our initial examination of the 
low lying spectrum appears to indicate the presence of light scalar bound states. 

We have also examined the issue of the phase of the Pfaffian induced by integrating over 
the fermions. Our Monte Carlo simulations have been conducted in the phase quenched 
ensemble in which this phase is neglected. We present evidence that the distribution of 
the phase measured in the phase quenched ensemble peaks around zero at large coupling. 
This corresponds to our expectation that the relevant Pfaffian is indeed real and positive 
definite in the continuum limit. Remarkably re- weighting with the phase in the usual way 
does not appear to generate large corrections to the phase quenched results at least in the 
case of the supersymmetric Ward identities. We have also examined the distribution of 
the eigenvalues of the scalar fields in the theory. The latter quantity yields information 
of the nature of the quantum moduli space. We provide strong evidence that the vacuum 
degeneracy of the classical theory is lifted at the quantum level. 

The paper is organized as follows; we first summarize the reformulation of the con- 
tinuum theory is terms of twisted fields. This twisted theory has a natural mapping to 
the lattice and we give the lattice action and the action of the twisted supersymmetry 
on the lattice fields. The details of our numerical algorithm are then described and lead 
to a presentation of our results both for U(2) and SU{2) gauge groups. A final section 
summarizes our conclusions and discusses future work. 

2. Continuum Twisted Formulation 

Consider a theory with AT = 2 supersymmetry in two dimensional Euclidean space. Such 
a theory contains two Majorana supercharges q l a transforming under the global symmetry 
group SO (2) x SO(2)r where the subscript corresponds to two dimensional rotations and 
the superscript describes the behavior under the R-symmetry corresponding to rotating 
the two Majorana fields into one another. 

The basic idea of twisting, which goes back to Witten 
new rotation group 

50(2)' = diagonal subgroup(50(2) x S0(2) R ) 

and to decompose all fields now as representations of this new rotation group. In simple 
terms what this means is that whenever I do a rotation in the base space by some angle I 
must do an equal rotation in the R-symmetry space. It is equivalent to treating the two 
indices I and a as equivalent. Thus the supercharges are to be interpreted as matrices in 
this twisted picture. 

qi -> Qap (2.1) 
Such a matrix may then be expanded on a basis of products of 2D gamma matrices 

q = QI + + Qi27i72 (2-2) 



27, 23, 28], is to introduce a 
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The fields (Q, Q^, Q12) are called the twisted supercharges. The original SUSY algebra then 
implies a corresponding twisted algebra which takes the form 

{^W=4 7 > M (2.3) 

In components this reads 

{Q,Q} = {Qi2,Qn} = {Q,Qu} = {Q„,Qu} = 
{Q,Q»} = Pfi 

{Qi2,Q»} = -e^Pv (2.4) 

Notice that the momentum is Q-exact and this fact suggests that generically both the 
entire energy- momentum tensor and hence the action will also be Q-exact. Furthermore, 
it should be clear the fermions admit a similar decomposition 

7] 

■qr = -I + ^ 7/1 + xi27i72 (2.5) 

and the twisted theory will not contain spinors but antisymmetric tensor fields. It is possible 
to abstract these p-form fields and consider them as components of a single geometrical 
object - a Kahler-Dirac field. 



* = (7^X12) (2-6) 



•1 
'2' 

A general complex Kahler-Dirac field describes two Dirac spinors in two dimensions. 
These spinors can be read off as the columns of the original fermion matrix in a particular 
7 matrix basis. However if we think of the twisted component fields as real then ^ = ^ T 
and the Kahler-Dirac field will describe two Majorana fields as required by the M = 2 
supersymmetry. Furthermore, any theory invariant under the scalar Q-symmetry must 
contain commuting superpartners for these fermionic fields. 

$ = $,A„,Bi 2 ) (2.7) 

The action for the Yang-Mills model may be written as S = f3Qk where 

A = Tr J d 2 x Qr?^, 0] + 2*12^12 + Xi2B 12 + (2.8) 

where each field is in the adjoint of some gauge group X = YIH^Lq 1 X a T a and we employ 
antihermitian generators T a . In practice we will consider both U(N) and SU(N) gauge 
groups with the convention that in the former case the traceless generators correspond to 
a = . . . N g — 2 and the final generator is proportional to the unit matrix T N " 1 
action of Q on the fields is given by 





= Yv 




= -Drf 


Q<i> 


= 


Qxn 


= Bu 


QB 12 


= y>,xi2] 




= V 


Qv 


= [<M 



(2.9) 
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Notice that Q 2 generates an (infinitessimal) gauge transformation parametrized by the field 
(f>. Carrying out the Q-variation and subsequently integrating out the multiplier field B\ 2 
we find the on-shell action 



S = (3Tr J d 2 x Q[<M 2 " \vl<f>,ri]-Fi2 



- D^D^cf) - Xi2[0>Xi2] 

+ (2-10) 

The terms involving the twisted fermion fields correspond to the component form of the 
Kahler-Dirac action 



y*.(d A -4)* ( 2 -n) 



where d>A is the gauged exterior derivative and ^ is the real Kahler-Dirac field introduced 
earlier (we will only consider flat space in this paper). 

3. Lattice Theory 

This continuum twisted action can be discretized according to the prescription detailed in 
Q. The on-shell lattice action is 

Sl = ^Y,(\[^] 2 + f1 2 F 12 

X ^ 

- -x + i2^xi2] (12) 

- 2^D+| + h.c) (3.1) 

where scalars such as (j)(x) are associated with lattice sites, vectors such as U^(x) with 
links and rank 2 tensors with plaquettes. For example, Xi2{%) is a lattice field associated 
with the 12-plaquette at site x. They are assigned corresponding gauge transformation 
properties [p6j] : 

(f>(x) -» G{x)(p{x)G- 1 {x) 
U^x) -> G(x)C/ M (x)G- 1 (x + ^) 

Xi 2 (x) -» G{x) X i2(x)G-\x + 1 + 2) (3.2) 

where G = is a lattice gauge transformation. Notice that for fields of non-zero spin 
the infinitessimal form of this lattice gauge transformation naturally leads not to the usual 
naive commutator characteristic of adjoint fields but to a shifted commutator. For link 
fields this looks like 

[ ( j>{x)f fX {x)-f lx {x)^x + ^} (3.3) 

In order to allow us to construct gauge invariant quantities using the above gauge transfor- 
mation rules it is essential that each continuum field with non-zero spin gives rise to two 
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lattice fields. This doubling of degrees of freedom can be associated to the two possible 
orientations of the underlying p-cube (for p > 0) on which the field lives. The doubling can 
be conveniently encompassed by promoting each field from real to complex and assigning 
the conjugate fields to transform in the following way 



Ul{x)^G(x + n)Ul{x)G-\x) 

X \ 2 (x) - G(x + 1 + 2) X t 12 (x)G- 1 (x) (3.4) 



This complexification has immediate consequences - the gauge group of the lattice theory 
is promoted from U(N) to GL(N,C) (or SU(N) to SL(N,C)) and the usual gauge links 
are no longer unitary matrices. In addition to the fields we need to define the lattice 



derivatives appearing in the lattice action eqn. 3.1. The action of the gauge covariant 



forward difference operator acting on scalars and vectors is defined by |2 

D+f(x) = U^x)f{x + M ) - f{x)U„{x) 

D+f u (x) = U^x)f v {x + /i) - U(x)U^x + v) (3.5) 

It clearly reduces to the usual gauge covariant derivative acting on adjoint fields in the 
naive continuum limit. Note that this derivative acting on a scalar or site field yields a 
field which gauge transforms as a link field and the corresponding derivative of a link field 
yields a field which transforms like a tensor or plaquette field. These properties allow us 
to construct a lattice analog of the continuum gauged exterior derivative by appropriately 
anti-symmetrizing in the spacetime indices. It is then possible to make a straightforward 
transcription of the continuum twisted action to the lattice. Furthermore, it is possible 
to write a covariant backward difference which is adjoint to the above operator for gauge 
invariant quantities. Its action on link and plaquette fields is given by 

D~U(x) = U(x)Ul(x) - Ul(x - fif^x - /i) 

D »Uv(%) = ffj, u (x)Ul(x + v)- L/J(z - n)fp,(x - n) (3.6) 



The action in eqn. |3.l| is invariant under the following lattice supersymmetry transformation 



which is a simple generalization of the continuum one given earlier 



QU^ = 




QVv = 




Q<t> = 





Qxn = 




QB 12 = 


[</>,Xl2] 


Q4> = 


V 


Qv = 





(12) 



(3.7) 

where the continuum D^cf) is replaced by the lattice forward difference D~^<p as required by 
gauge invariance and the prime on the commutators reflects its shifted nature as discussed 
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earlier. The transformations of the conjugate fields are gotten by taking the adjoint of 
these variations with the constraint that f< = —f for any scalar or site field. Notice that it 
is only for the group GL(N, C) that the adjoint of a link field can be taken as transforming 
independently from the link field itself. Finally the Yang-Mills field strength appearing 
above is defined by 

Ffiv (x) = D+U v (x) -» F™ nt as a - (3.8) 
Writing this term out we find 

0T I J2f! 2 (x)F 12 (x) (3.9) 

X 

pTr ( 2I -Up- U P \ + /3Tr £ (Af 12 + M 21 - 21) (3.10) 

X X 

where 

Up = U x {x)U 2 {x + l)u\{x + 2)ul{x) (3.11) 
resembles the usual Wilson plaquette term and 

M 12 = U^x)ul{x)ul(x + l)U 2 (x + 1) (3.12) 

is a new zero area Wilson loop term which would vanish if the link variables were restricted 
to unitary matrices. Notice the appearance of the Wilson term depends crucially on the 
appearance of both and F^ v which lends some support to the use of complex variables 
in the formulation of the theory. 

The lattice action we have written down possesses one additional SO(l, 1) symmetry 
corresponding to the transformations 

ipn^Xipn, V, X12 -» \r}, xi2 (3.13) 



1 - 



\ 2 <p (3.14) 



The transformation of the conjugate fields under this global symmetry is identical. This 
symmetry is useful as it guarantees the absence of additive mass renormalizations in the 
lattice theory. 

Finally we should point out that the spectrum of this lattice theory contains no lattice 
doubles either fermionic or bosonic. This result follows from the work of Rabin, Becher 



and Joos 30, 31] who show that actions written in terms of exterior derivatives and 
tensor fields may be discretized without generating doubled modes. The discretization 
prescription is given explicitly by replacing the usual partial derivatives in the continuum 
theory by appropriate difference operators: 

— > D + if acts like d 
—fD~ii acts like S 



-7- 



This double free property can be seen explicitly in our case by examining the form of the 
fermion action. Using the following decomposition of the Kahler-Dirac field 



* = 



Xl2 



(3.15) 



■01 
\4>2 ) 



Our twisted fermion action can be recast in the form ^ M(U,4>)^ where the fermion 
operator M is given in block form by 



with the Yukawas lying along the diagonal and K the lattice Kahler-Dirac operator taking 
the form 



In the continuum theory the Kahler-Dirac field satisfies a reality condition ^ = ^> T and 
integration over these anticommuting fields yields the Pfaffian of the fermion operator 
Pf(M). In the free limit we find that this prescription, when applied to the above lattice 
operator, yields the determinant of an explicitly double free lattice laplacian Pf(M) = 



As an aside we note that there is well known equivalence between Kahler-Dirac fermions 
and staggered fermions - the 4 component fields of a single Kahler-Dirac field in two 
dimensions can be mapped to site fields on a lattice of half the lattice spacing and the 
free Kahler-Dirac action goes over into the usual staggered action. This is another way 
of understanding why discretizations of the Kahler-Dirac theory avoid spectrum doubling. 
The usual flavor replication of staggered fermions here becomes a bonus - it yields an 
automatic description of the two degenerate fermions required by N = 2 supersymmetry. 
Of course it must be remembered that the gauging of our lattice Kahler-Dirac action is not 
at all the usual gauging of staggered fermions so the exact equivalence does not persist in 
the interacting theory. 

To conclude this description of the lattice theory we should return to the issue of 
complexification. The lattice formulation we have given requires a doubling of degrees of 
freedom - we have argued that this is quite natural in any lattice theory and can be asso- 
ciated with the two possible orientations of the underlying p-cube. We have parametrized 
this doubling in terms of complex fields. However, the target continuum theory that we 
are hoping to reproduce in the limit of vanishing lattice spacing corresponds to putting the 
imaginary parts of the fields to zero or more accurately to setting 




(3.16) 




(3.17) 



det(tf) = det(D+D~) 




(3.18) 
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Let us examine what this means for the complexified lattice theory. First, consider the 
effective action that results from integrating out the grassmann variables. While the com- 
plexified theory would lead to a determinant, integration in the truncated theory should 
result in a Pfaffian. If we ignore possible phase problems and replace the latter by a square 
root of the determinant we can see that the fermion effective action will still be gauge 
invariant after truncation to the real line. The bosonic action also remains gauge invariant 
after the projection and clearly both contributions target the correct continuum theory in 
the classical continuum limit. The remaining important issue is whether the Ward iden- 
tities corresponding to the twisted supersymmetry still hold in the truncated theory (or 
more conservatively, hold in the limit of vanishing lattice spacing with no additional fine 
tuning). We conjecture that this may be so and have followed this approach so far in our 
numerical work. 

It is possible to make some progress in understanding why this might indeed be true. 
First parametrize the general GL(N,C) link field U^ix) in terms of a unitary component 
u ll (x) and a positive definite hermitian piece R^(x) in the following way 

U^x) = R^x)u^x) (3.19) 

Now consider the second term in the gauge action 

/3Tr ( M i2 + M 2 i - 21) (3.20) 

X 

and insert the general decomposition of the gauge link given above. The result for M12 is 

M12 = R\{x)Rl(x + l) (3.21) 

with a similar result for M21. Notice it is independent of the unitary piece u^{x). Consider 
the theory in the continuum limit [3 — > 00. It should be clear that in such a limit each R^x) 
is driven to the identity and the complex bosonic action smoothly approaches the usual one 
involving real fields. Furthermore, for large f3 the fermion operator is both independent of 
Rfj,(x) and antisymmetric. Thus in this limit the real and imaginary components of the 
fermions decouple and the fermion determinant factors into the square of a Pfaffian. This 
decoupling ensures that expectation values of operators depending only on the real part of 
the fermion field (the Majorana condition) and computed for large (5 will approach their 
values computed in the truncated theory. Thus, the Ward identities of the truncated theory 
should hold at least for large (3 since they can be viewed as coming from the complexified 
theory (possessing explicit exact Q-symmetry) in the limit of infinite (3. Of course large 
(3 also corresponds to the limit of vanishing lattice spacing and we see that this line of 
reasoning constitutes an argument that the Ward identities in the truncated lattice theory 
should be realized without fine tuning in the continuum limit. As we shall show in the 
next section our numerical results are consistent with this. Indeed, in the case of SU (2) 
the Ward identities appear to hold with small errors even for small f3. 
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4. Simulations 



For our simulations we have taken the gauge links to be unitary matrices taking values 
in either £7(2) or SU(2). As in the continuum theory the scalars <f> and (j) are taken to 
be complex conjugates of each other 3 . In the case of £7(2) the bosonic action possesses an 
exact zero mode 4>o = (0, 0, 0, cf) 3 ) which we hence regulate with the addition of a mass term 
m 2 (f>(x)(j)(x). At the end of the calculation we should like to send m — ► to recover 
the correct target theory. The bosonic action is real positive semi-definite, gauge invariant 
and clearly has the correct naive continuum limit. 

To this gauge and scalar action we should add the effective action gotten by integrating 
over the grassmann valued fields. We will represent this as 

det3(M(£7» t M(C/,0) + m 2 ) (4.1) 

where M is the lattice Kahler-Dirac operator introduced earlier. The power of \ reflects 
the Majorana nature of the continuum Kahler-Dirac field. Notice that we have added a 
gluino mass term for the fermions which regulates the corresponding fermion zero mode 

= (0,0, 0, ry 3 ) arising in the £7(2) theory as a consequence of supersymmetry. The 
SO(l, 1) symmetry in the massless case prohibits additive renormalization of this mass as 
a result of quantum effects. In the case of SU{2) this mass parameter can be set to zero. 
Clearly, the form of the fermion effective action we employ does not take into account 
any nontrivial phase associated with the fermion determinant or Pfaffian - our simulations 
generate the phase quenched ensemble. We later examine the phase explicitly. 

To simulate this model we have used the RHMC algorithm developed by Clark and 
Kennedy [32]. The first step of this algorithm replaces the effective action by an integration 



over auxiliary commuting pseudofermion fields F, in the following way 

deti(MtM + m 2 ) = J DFDF^ e~ F ^ M " M+m2 )~ k F (4.2) 

The key idea of RHMC is to use an optimal (in the minimax sense) rational approximation 
to this inverse fractional power. 

1 P{x) 



xi Q( x ) 

where 



(4.3) 



N-l N-l 



P(x) = £ Q(z) = E * x< ( 4 - 4 ) 

8=0 8=0 

Notice that we restrict ourselves to equal order polynomials in numerator and denomi- 
nator. In practice it is important to use a partial fraction representation of this rational 
approximation 

1 N 

X4 i=1 X T Pi 



3 the antihermitian nature of our basis T a actually ensures that <j>' = 
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The coefficients Oj, for i = 1 . . . N can be computed offline using the remez algorithm 4 . 
Furthermore, the coefficients can be shown to be real positive. Thus the linear systems 
are well behaved and unlike the case of polynomial approximation the rational fraction ap- 
proximations are robust, stable and converge rapidly with N. The resulting pseudofermion 
action becomes 

N 

SpF=aoF t F+ |>t __^_ f («> 

It is thus just a sum of standard 2 flavor pseudofermion actions with varying amplitudes and 
mass parameters. In principle this pseudofermion action can now be used in a conventional 
HMC algorithm to yield an exact simulation of the original effective action fl33[| . This 
algorithm requires that we compute the pseudofermion forces. For example, the additional 
force on the gauge links due to the pseudofermions takes the form 

i=l 

where the vector \ l is the solution of the linear problem 

+ m 2 + = F (4.8) 

The final trick needed to render this approach feasible is to utilize a multi-mass solver to 
solve all N sparse linear systems simultaneously and with a computational cost determined 
primarily by the smallest shift We use a multi-mass version of the usual conjugate 
gradient CG algorithm |j34| . In practice for the simulations shown here we have used 
N = 15 and an approximation that gives an absolute bound on the relative error of 10 -6 
for eigenvalues of M'M ranging from 10 -8 to 10 which conservatively covers the range need 
for both our U(2) and SU (2) runs. We monitor the spectrum continuously to make sure 
our approximation remains good. Typically we have amassed between 10 3 and 10 4 HMC 
trajectories for each set of parameter values which leads to statistical errors of between 0.1 
and 2.0 percent depending on observable. 

Finally we make some remarks on the representation of this fermion operator. For 
the purposes of computation this abstract lattice fermion operator M is replaced with a 
sparse matrix whose non-zero elements are gotten by choosing an explicit basis for the 
group generators and evaluating all traces over internal indices. For example the term 
Y.x Tl ^U x ) D ^v(x) yields 

^(x)V^(x)r,\x + /i) - ^(x)V^(x)r, b (x) (4.9) 

where V^ b (x) = Ti(T a V " /i (x)T b ) . Similarly Yukawa terms reduce to matrix elements of the 
form 

rf(x)(<S> ab (x) - $ ba (x))ri\x) (4.10) 
where <J> a6 (x) = Tr(T a <f)(x) c T c ) T b ). In practice we use the basis 

T° = iai T 1 = ia 2 T 2 = ia 3 (4.11) 
4 many thanks to Mike Clark for providing us with a copy of his remez code 
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where <Tj are the usual Pauli matrices and in the case of U(2) the generator T 3 = il 

As usual the running time of the simulations is dominated by the need to solve the 
linear system eqn. [D^ for each step down a Monte Carlo trajectory. We use the usual 
sparse matrix techniques to optimize the CG-solver to accomplish this. 

5. Phase quenched U(2) model 

The numerical results we present here come from simulations where the lattice length takes 
the values L = 2, 4, 8 for mass parameter m = 0.1 and L = 3, 6 for m = 0.01. The coupling 
varies over the range (3 = 0.5 — ► 4.0. Clearly Ward identities corresponding to the twisted 
supersymmetry Q are of prime interest. They are simply expectation values of the form 
< QO > and should be zero by supersymmetry. Perhaps the simplest of these corresponds 
to the action itself < S = QA >=< 5b > + < Sf >= 0. This fact, together with the 
quadratic nature of the fermion action allows us to compute the bosonic (gauge plus scalar) 
action exactly using a simple scaling argument and for all values of the coupling constant 
j3 we find 

(3 < S B >= ^N g L 2 (5.1) 

where N g = 4 is the number of generators of U{2). Figure 1. shows a plot of < Sb > for 
a range of couplings (3 and three different lattice sizes using a mass m = 0.1. The bold lines 
shows the analytic prediction based on supersymmetry (for clarity we have added multiples 
of —0.25 to the curves and lines to split up the data from different lattice sizes) While there 
are clearly deviations of order 4-5% at small coupling these appear to disappear at large (3 
in line with our expectations. Figure 2. shows equivalent data for L = 3, 6 at the smaller 
gluino mass m = 0.01. Again the horizontal lines show the analytic result expected from 
supersymmetry. In this case the deviations for the small lattice appear larger at small 
coupling but nevertheless appear to converge toward the theoretical expectation on the 
basis of supersymmetry as (3 is increased. The larger lattice L = 6 data are even better. 
We have also examined other Ward identities corresponding to the local operator choices 

O = Ox{x) =rj(x)[(f>(x),0(x)} 
O = 2 {x) = X^F^ 

= 3 (x)=^^(x)D+^(x) (5.2) 

After Q-variation we find 

QOt = [<f>,W-v[<t>,v] 

Q0 2 = Fl v {x)F^{x) -x^D+ij} v] 

QOs = - ^D+ri - VvfVv, 4>] (5.3) 

The results for expectation values of these Q-variations for (3 = 4.0 and m = 0.01 are 
displayed in tables 1. and 2. corresponding to lattice sizes L = 3 and L = 6 respectively 
(we denote the bosonic contribution to the Ward identity by B and the fermionic one 
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o 




Re(F) 


Im{F) 


Oi 


0.88(7) 


-0.86(6) 


-0.0003(5) 


o 2 


0.285(16) 


-0.297(16) 


-0.0009(28) 


o 3 


0.90(45) 


-0.986(4) 


0.0035(39) 



Table 1: Ward identities for U{2) and (3 = 4.0,L = 3 






B 


Re(F) 


Im{F) 


Oi 


0.38(8) 


-0.43(10) 


-0.0008(15) 


o 2 


0.409(22) 


-0.386(25) 


-0.0006(7) 


03 


1.22(49) 


-0.995(1) 


-0.0003(20) 



Table 2: Ward identities for U{2) and (3 = 4.0,L = 6 



by F). Notice that the imaginary part of the fermionic correlator is always small and 
statistically consistent with zero 5 . Within the statistical errors the bosonic and fermionic 
contributions do add to zero confirming the presence of the Q-symmetry in the quantum 
lattice theory for this coupling. 

The apparent breaking of supersymmetry at small (3 appears to be correlated to the 
symmetry properties of the lattice fermion operator M(U,(j)). In the continuum this op- 
erator is complex antisymmetric and hence its determinant can be written as the square 
of a Pfaffian. This is true for both SU(2) and U(2) theories (the U(l) or trace mode of 
the scalars and gauge field disappear from the fermion operator for all couplings rendering 
both continuum theories equivalent in this respect). 

The antisymmetric condition is very important also in the lattice theory because in 
this case we have argued that the real and imaginary components of the fermion decouple 
- a necessary condition for existence of a Pfaffian and the truncation of the theory to the 
real line. If the matrix is not antisymmetric this factorization cannot be achieved and 
supersymmetry must necessarily be broken. In the case of SU{2) it is not hard to see that 
the matrix is indeed antisymmetric and we exploit this fact later when we examine the 
phase of the Pfaffian for the SU{2) system. This property is not shared by the U{2) model 
in general though since the trace degrees of freedom of the gauge link do not decouple for 
finite (3. Approximate decoupling will now occur in the region of large coupling (3 where 
the approximation U^{x) = 1 + A^(x) + • • • becomes accurate. By examining the fourth 
root of the plaquette we learn that \A^\ ~ 0.05 for (3 = 3.0. Hence corrections to this 
approximation will be less than one percent if (3 > 3 consistent with the restoration of 
supersymmetry we see in the large (3 regime. This decoupling can be seen explicitly in 
figure 3. which shows the Monte Carlo evolution of the scalar field for (3 = 3.0, L = 6 
and m = 0.01 in the U(2) theory. Two modes are shown corresponding to (j) 3 (the trace 

5 the fermionic correlator < >= \M~^ where spacetime, group and Kahler-Dirac indices are 

combined into a single index i and the factor of | originates from the Majorana nature of the Kahler-Dirac 
field 
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mode) and (fP a traceless mode. The trace mode <fi 3 behaves as a quasi massless degree of 
freedom undergoing large fluctuations regulated only by the imposed IR cut-off of m = 0.01 
while the traceless degrees of freedom fluctuate independently over a scale two orders of 
magnitude smaller. The spectrum of a typical equilibrated configuration for j3 = 4.0 on a 
lattice of size L = 6 is shown in figure 4. Much of the spectrum is concentrated close to 
the imaginary axis and is indicative of light continuum-like states. Notice the approximate 
pairing of eigenvalues (A, —A) related to the approximate antisymmetry of the fermion 
operator at this large coupling. There are in addition "islands" of additional states at 
large eigenvalue which we conjecture are related to the large fluctuations of the nonzero 
trace modes of the scalars. These don't appear to have any continuum interpretation. 

6. Phase quenched SU(2) model 

Since the SU (2) theory contains no exact zero modes we have been able to simulate the 
model at exactly zero gluino mass. Figure 5. shows a plot of the mean bosonic action 
normalized to unity as a function of (3 = 0.5 — > 4.0 for L = 3, 6, 8 for the SU (2) theory. In 
contrast with the £7(2) lattice theory the scalar supersymmetry appears to be good here 
down to small coupling (3. As we have remarked we conjecture that this is related to the 
presence of an exact antisymmetry of M in the SU (2) case. Clearly the Yukawas possess 
this property. What is non trivial is that this is also true of the gauged lattice Kahler-Dirac 
term. In general this term is antihermitian but in the case of SU{2) it is also real. This 
follows from the special property of the Pauli matrices 

e^ a ' = cos W + ^J^U<7i (6.1) 

H 

Using this representation it is easy to show that V® b (x) is real. Figure 6. shows a plot of 
the eigenvalues of the SU{2) theory for L = 6 at both (3 = 0.5 and (3 = 4.0 in which the 
exact pairing of eigenvalues is manifest. Notice that as (3 increases the real parts of these 
eigenvalues decrease. In the limit we expect the eigenvalues to lie along the imaginary 
axis yielding a real, positive definite determinant as for the continuum theory. We have 
additionally measured the same local Ward identities as for J7(2) with the results listed in 
tables 3. and 4. Here, the data is taken from runs with L = 6 and two different values 
of the coupling (3 = 0.5 and (3 = 4.0. Again for both small and large coupling these local 



o 


B 


Re(F) 


Im(F) 


O 


B 


Re(F) 


Im(F) 


Oi 


4.70(20) 


-4.72(04) 


-0.04(04) 


Oi 


0.167(10) 


-0.195(8) 


-0.007(07) 


o 2 


1.95(07) 


-2.02(02) 


0.014(14) 


o 2 


0.344(10) 


-0.346(1) 


-0.0007(8) 


o 3 


5.76(15) 


-5.74(05) 


-0.043(35) 


o 3 


0.759(21) 


-0.739(4) 


0.0012(17) 



Table 3: Ward identities for SU(2) and (5 = Table 4: Ward identities for SU(2) and (3 = 
0.5,L = 6 4.0,L = 6 



Ward identities appear to be satisfied to within statistical error. 

Since the simulations of the SU (2) are carried out at zero mass we have a priori no lower 
bound on the eigenvalue spectrum and so we monitor the smallest eigenvalue continuously 
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to ensure it lies within the limits required by our minimax approximation to the inverse 
fourth root of the fermion operator 6 . A typical plot of the Monte Carlo evolution of |A m i n | 
is shown in figure 7. We observe that the magnitude of this smallest eigenvalue decreases 
with increasing (3 and L but we see no evidence for an exact zero mode as would be expected 
for a supersymmetric theory whose flat directions survive quantum corrections. We will 
return to this issue when we discuss the quantum moduli space. 

Since we employ periodic boundary conditions the partition function we simulate yields 
the Witten index of the theory and is explicitly independent of coupling constant (recall 
that ^0- =< 5 >= 0). This in turn implies that there can be no thermodynamic singu- 
larity for finite (3 and we expect the theory exists in a single phase. Figure 8. confirms 
this expectation by plotting the string tension as estimated from the 2x2 Creutz ratio as 
a function of (3 for a lattice of size L = 6. The string tension appears to be non-zero and 
smoothly varying over this range of coupling. 

We have also examined a couple of correlation functions which give us direct access to 
the low-lying mass states of the theory. The simplest 50(1, 1) and gauge invariant bosonic 
correlator takes the form 

G B (t) = Y,< ^ a ( x > t)T(x, t)cj) b {x', 0$V, 0) > (6.2) 

x,x' 

and the sum over spatial sites x, x' projects to the zero momentum sector. We have also 
examined a fermionic correlator of the form 

G F (t) = ^2< rf{x, t)(f> a (x, t)rj\x', 0)O V, 0) > (6.3) 

x,x' 

where O = [cf>, </>]. These functions are shown in figure 9. together with fits to hyperbolic 
cosines. As a consequence of supersymmetry we expect that the lowest lying bosonic and 
fermionic states should have the same mass. Within statistical errors this is true (and 
indeed this state is rather light). However the errors on the fermion are large and the 
current data is really inadequate to decide this question. Of course it is also not clear we 
have the correct interpolating operator for the lightest fermion state - further investigations 
of these issues are underway. 

7. Reality of the fermionic effective action 

Up to this point we have neglected a possible phase associated with the Pfaffian induced 
by integration over the fermion fields. As we have seen the truncation to the real line 
in general breaks the supersymmetry in the U (2) model so we will concentrate on the 
SU(2) case where the fermion operator is an antisymmetric matrix and a Pfaffian can be 
unambiguously defined. As usual we can always compensate for neglecting this phase in 
the Monte Carlo simulation by re-weighting all observables by the phase factor according 
to the simple rule 

< o >= ' " =0 (7.1) 

< e ia ( u '<t>) > a=0 



'For the U(2) theory this is rigorously bounded below by m 
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We thus have computed the Pfaffian as one of our observables allowing us to carry out this 
re-weighting procedure when computing expectation values. The Pfaffian computation is 
carried out by using a variant of Gaussian elimination with full pivoting to transform the 
2n x 2n dimensional antisymmetric matrix M into the canonical form 



(7.2) 
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Then 



Pf (M) = H Xi 



(7.3) 



Consider first the phase itself. Figure 10. shows a plot of the distribution of cos a for 
L = 3 and (3 = 4.0. A strong peak close to a = is manifest. This peak strengthens with 
increasing (3 as the scalars (which are responsible for a non-zero phase) are driven closer 
to zero. Figure 11. shows a plot of cos a vs (3 for L = 6 (< sin a > is small and always 
statistically consistent with zero). We see that it increases from values close to zero to 
attain 

<cosa>~0.2 /3 = 4.0 (7.4) 

Hence, in the range of coupling we have simulated, it clearly fluctuates strongly from the 
naive value of unity used in generating the phase quenched ensemble. In light of this we 
have re-examined the Ward identities now weighted with this phase factor. Table 5. shows 
the mean re-weighted bosonic action together with the phase quenched value for L = 3 and 
all (3 (note that these numbers are not normalized to unity as in the earlier plots). While 
re-weighting typically amplifies the estimated error it does not appear to change the mean 
value for this observable at least within the statistical errors. 

This conclusion is strengthened by examin- 
ing other re-weighted Ward identities correspond- 
ing to the set of operators Oi, O2 and O3 given 
earlier. Tables 6. and 7. compare the naive 
(phase quenched) expectation values with their re- 
weighted values for lattice size L = 6 and [3 = 0.5 
and (3 = 4.0. Again, there is no evidence that 
the central values change within the (admittedly 
large) statistical errors. 

Taken at face value this apparent weak de- 
pendence of the expectation values on reweight- 
ing seems to indicate that the phase fluctuates ap- 
proximately independently of the other observables leading to an, at least approximate, 
factorization in the reweighted observable 



p 


S B 


Reweighted Sb 


0.5 


40.35(19) 


39.80(95) 


1.0 


40.45(18) 


40.42(117) 


2.0 


40.11(32) 


40.80(130) 


2.5 


40.16(32) 


40.62(130) 


3.0 


39.74(25) 


40.11(140) 


3.5 


39.99(30) 


40.58(200) 


4.0 


39.73(30) 


39.66(220) 



Table 5: Ward identities for SU(2) and 
P = 4.0,L = 6 



< O > fuU = 



< Oe ia > a=0 <0> c 



=o< e > a =0 



< e ia >a=0 



< e ia >a=0 



=< o > a= 



(7.5) 
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o 


B 


Re{F) 


Reweighted B 


Reweighted Re(F) 




4.70(20) 


-4.72(03) 


6.46(210) 


-5.04(100) 




1.95(07) 


-2.02(02) 


2.14(100) 


-2.01(50) 


o 3 


5.76(15) 


-5.74(02) 


5.74(200) 


-5.25(200) 



Table 6: Naive vs reweighted Ward identities for SU(2) and (5 — 0.5,i = 6 






B 


Re(F) 


Reweighted B 


Reweighted Re(F) 


Oi 


0.167(10) 


-0.195(7) 


0.159(40) 


-0.137(50) 


o 2 


0.344(10) 


-0.346(1) 


0.356(60) 


-0.348(50) 


o 3 


0.759(21) 


-0.739(4) 


0.717(110) 


-0.733(120) 



Table 7: Naive vs reweighted Ward identities for SU(2) and j3 = 4.0, L = 6 



As a practical matter this means that expectation values obtained within the phase quenched 
approximation may be quite reliable in spite of the large phase fluctuations. 

8. Quantum moduli space 

Finally we turn to an important issue concerning the two scalar fields that appear in this 
theory. The classical vacua allow for any set of scalars which are constant over the lattice 
and satisfy 

[<M=0 (8.1) 

In the case of SU (2) and using the parameterization <\> = fa + ifa we find vacuum states 
of the form 

fa = (A,0,0) fa = (B,0,0) (8.2) 

together with global SU{2) rotations of this configuration. Thus we have a classical vacuum 
state for any value of A and B. This space of vacuum solutions is referred to as the moduli 
space of the theory. The presence of such a non-trivial moduli space corresponds to the 
existence of flat directions in the theory. This is problematic in the quantum theory as 
integration over such flat directions may induce IR divergences. However we find that 
this is not the case in practice - the quantum ground state appears to be unique and the 
flat directions are lifted (except for the trivial U(l) factor associated with trace part of 
U(2)). This can be seen from the distribution of eigenvalues of the scalars averaged over 
the lattice. Figure 12. shows a plot of this distribution for the SU(2) theory on a L = 3 
lattice both for (3 = 0.5 and (3 = 4.0. The distribution is symmetric about the origin so 
we show only the positive values here. Both for small and large coupling the distribution 
possesses a well defined peak with a tail extending out to large eigenvalue. Notice that the 
peak moves to smaller values as f3 increases in line with the observed suppression of scalar 
field fluctuations with increasing coupling. The data indicates that the partition function 
at least exists and most likely at least some of the lower moments of the scalar field. 
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This result is reminiscent of similar results obtained for zero dimensional supersymmetric 
Yang- Mills integrals [35]. 



9. Conclusions 

In this paper we have presented initial results from a full simulation of the M = 2 super 
Yang-Mills theory in two dimensions. The lattice action we employ was derived in Q and 
follows from a reformulation of the theory in terms of twisted fields. It is invariant under 
a global 50(1, 1) symmetry and U(N) lattice gauge transformations and exhibits, at least 
in its complex form, an exact scalar supersymmetry. We show results for both U{2) and 
SU(2) theories for a range of lattice size L and coupling (3. To check for supersymmetry 
we have examined a number of Ward identities. 

In the case of U{2) the supersymmetry appears to be only exact at large (3. At 
small (3 the fermion operator is not antisymmetric, its Pfaffian is not even defined and 
the truncation of the complexified theory to the real line appears to break supersymmetry. 
However, the theory does not seem to need any fine tuning to regain supersymmetry for 
large coupling (3 and hence in the continuum limit. 

In the case of SU{2) the lattice gauged Kahler-Dirac operator is real and antisymmetric 
and supersymmetry is manifest for all couplings - we see no statistically significant violation 
of any Ward identities at the 1% level for any coupling or lattice size. We conjecture that 
the key property which allows the supersymmetry to be realized is the antisymmetric 
property of the fermion operator. Since the latter is always antihermitian this in turn boils 
down to a reality property on the gauged Kahler-Dirac operator. At first glance it appears 
that SU{2) is rather special in this regard. However, reality of the SU(N) theories will be 
guaranteed for all N if instead of using the N x N basis for the generators we employ the 
basis constructed from the structure constants themselves Tg c = f a t, c . Such a choice yields 
the same naive continuum limit but guarantees that the gauged Kahler-Dirac action is real 
and antisymmetric. We conjecture that such models will resemble SU(2) and exhibit an 
exact scalar supersymmetry at the quantum level. 

In the case of SU{2) we have also examined the validity of the phase quenched ap- 
proximation used in our simulations. The phase appears to approach zero in the large j3 
continuum limit as naively expected. Remarkably, re-weighting observables with the phase 
does not appear to have a strong influence on expectation values even for small f3 at least 
in the case of observables corresponding to Ward identities. We also show results on the 
distribution of the eigenvalues of the scalar fields. This distribution possess a well-defined 
peak which narrows and moves to smaller values as (3 increases. This structure is similar 
to the case of zero dimensional SUSY Yang- Mills and indicate that the classical vacua are 
lifted via quantum effects. 

To summarize, our initial numerical investigations of twisted formulations of lattice 
super Yang-Mills theories are quite positive - it appears that full dynamical simulations of 
these theories are possible with rather moderate computational resources (the work pre- 
sented here was obtained with 0(300) single CPU days on a P3 cluster). These preliminary 
results provide evidence that supersymmetry is indeed realized at the quantum level and 
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have already allowed us to check certain qualitative features of the theory - for example, 
our data support the lifting of the classical vacua and finiteness of the partition func- 
tion. They also support a single phase picture for this theory. It would be nice to extend 
these calculations to larger lattices to be able to get more reliable results for the low-lying 
spectrum which could be compared to the SDLCQ results reported in [36] and to the di- 
mensionally reduced quantum mechanics case |37j. It is important also to examine Ward 
identities corresponding to other elements of the twisted supersymmetry to see whether 
indeed the other super symmetries are realized without fine tuning in the continuum limit. 
Results from these investigations will be published elsewhere [38]. 

It should be stressed that low dimensional super Yang-Mills theories are of great in- 
terest because of their conjectured connections to various types of (super)gravity theory. 
They are a place where lattice simulations could potentially play an important role since 
high precision exact dynamical simulations are possible on large lattices. In principle, the 
strongest connections to gravitational systems are exhibited for theories with sixteen su- 
percharges rather than the four supercharge case considered here. However, dimensional 
reduction of the Af = D = 4 lattice action constructed in [24] would yield Q-exact actions 
for these systems whose fermion content could be represented using Kahler-Dirac fields. 
Notice that although the Kahler-Dirac action is related to the usual staggered fermion 
action there is no "fourth root" problem with these theories since the fermion degeneracy 
associated with these lattice actions precisely accounts for the number of physical fermions 
required by the extended supersymmetry. 
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Figure 1: gf^ < 5 B > for U{2) at to = 0.1 
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U(2), m=0.01 




Figure 2: < S B > for U{2) at to = 0.01 
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Figure 3: Monte Carlo evolution of trace and traceless components of scalars 
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Figure 4: Eigenvalue spectrum for {/ (2) theory at /3 = 4.0 
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Figure 5: M, < S B > for SU(2) 
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Figure 6: Eigenvalue spectrum for SU (2) theory at L = 6 
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Figure 9: Bosonic and fermionic correlators in SU{2) for L = 6 and /3 = 3.0 
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Figure 10: P(cosa) vs cos a for L = 3 [3 = 4.0 
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